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We analyze the dynamics of a Dirac-Born-Infeld (DBI) field in a cosmological set-up which includes 
a perfect fluid. Introducing convenient dynamical variables, we show the evolution equations form 
an autonomous system when the potential and the brane tension of the DBI field are arbitrary 
power-law or exponential functions of the DBI field. In particular we find scaling solutions can exist 
when powers of the field in the potential and warp-factor satisfy specific relations. A new class 
of fixed-point solutions are obtained corresponding to points which initially appear singular in the 
evolution equations, but on closer inspection are actually well defined. In all cases, we perform 
a phase-space analysis and obtain the late-time attractor structure of the system. Of particular 
note when considering cosmological perturbations in DBI inflation is a fixed-point solution where 
the Lorentz factor is a finite large constant and the equation of state parameter of the DBI field is 
W = — 1. Since in this case the speed of sound c s becomes constant, the solution can be thought to 
serve as a good background to perturb about. 

PACS numbers: pacs: 98.80.Cq 



I. INTRODUCTION 

The inflationary paradigm remains to date the most 
successful explanation for the origin of the observed 
temperature fluctuations of the cosmic microwave back- 
ground (CMB) (see, e.g., [H, for reviews). However, es- 
tablishing its origin in fundamental theory has not been 
quite as successful and so for many it remains a fasci- 
nating paradigm in search of an underlying theory. The 
favourite candidate for this is String theory and for a nice 
review of the construction of inflation models in string 
theory, see 

One interesting model, recently proposed from string 
theory is Dirac-Born-Infeld (DBI) inflation where 
inflation is driven by the motion of a D3-brane in a 
warped throat region of a compact internal space. In 
this model, since the inflaton is the position of a D-brane 
with a DBI action, its kinetic term is inevitably non- 
canonical. In addition to this kinetic term, its effective 
action includes a potential arising from the quantum in- 
teraction between D-branes, with the brane tension en- 
coding geometrical information about the throat region 
of the compact space. Because of these novel ingredients, 
the predictions of DBI inflation are quite different from 
the ones from the standard slow-roll inflation models and 
it has led to an intense period of research into the sce- 
nario, including work into the background dynamics and 
linear perturbations [9l4lfj|. 

From the phenomenological viewpoint, one of the main 
reasons why this model has attracted attention is be- 
cause of its sizable equilateral type of primordial non- 



Gaussinaity, first pointed by [4| and explored further in 
0-0, Il7l - t34l | . Such a large signal can not be obtained in 
standard slow-roll inflation, hence it opens up the possi- 
bility of distinguishing this model from that of the stan- 
dard slow-roll inflation- although not necessarily distin- 
guishing it from more general slow-roll models going be- 
yond single field inflation as they may also lead to large 
signals. Furthermore, because of the dependence of the 
effective four dimensional string tension on the nature of 
the compact internal space, it is possible to place further 
constraints on the parameters related with compactifica- 
tions such as flux numbers. For the most up to date ob- 
servational constraints on and consequences of DBI infla- 
tion, see [35M48I ] . Recently the idea of low scale inflation 
arising from the DBI action has been invoked to explain 
the late time acceleration we associate with dark energy 
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Given the results mentioned above and inspired by the 
possibility of inflation being an attractor solution in DBI 
models, we believe there is a need to understand the late- 
time attractor structure for as general a DBI set-up as 
possible. It is well known that for a canonical scalar field 
with a potential, scaling solutions can exist where the ra- 
tio of the kinetic and potential energy of the scalar field 
maintain the same ratio, and understanding the stability 
of these solutions is important in determining the na- 
ture of the late-time solutions [53l - l64l |. Among the ear- 
lier work, the phase-space analysis proposed by (57j is 
particularly powerful because it allows us to make use 
of suitable dimensionless dynamical variables, in order 
to establish the global stability of such scaling solutions. 
For related work which analyzes the dynamics in DBI 
models with general inflationary potentials, see [65l - [68j . 

Now, this method for analysing the stability of the so- 
lutions has only been applied to the case where the DBI 
field has a quadratic potential and the associated D-brane 
is in the anti-de Sitter throat 1691 . The late-time attrac- 
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tor structure for the case without a potential, known as 
tachyon cosmology, has been studied extensively, [70l - [78| 
and the case with another degree of freedom, has been 
studied in [zl|8(|. 

In this paper, in order to make the cosmological appli- 
cation of scaling solutions in DBI models more complete 
and as a natural extension of [69J , we obtain the late-time 
attractor structure of the system including a perfect fluid 
plus a DBI field whose potential and brane tension are 
arbitrary power-law or exponential functions of the DBI 
field. 

The rest of the paper is arranged as follows. In section 
IrEI we present the model and basic equations. Then, in 
section Mil we consider the models where the potential 
and brane tension are arbitrary power-law functions of 
the DBI field. Special emphasis is given to a new set of 
fixed-point solutions which at first site appear singular in 
the equations of motion, but on closer inspection are well 
defined. This is followed in section IIVI with an analysis 
of the models where the potential and brane tension are 
arbitrary exponential functions of the DBI field. Finally, 
we summarise in section Ivl 



II. BASIC EQUATIONS 

We consider a DBI field, <j), with the following effective 
action Q: 

s = - I dW^-^Vi - 2/(p - 1) - V(fi), 



where 



X 



(1) 



(2) 



V{4>) is a potential that arises from quantum interac- 
tions beween a D3-brane associated with <f> and other D- 
branes. Although a quadratic potential was considered 
in discussions on the exact form of the potential are 
still ongoing. f(4>) is the inverse of the D3-brane tension 
and contains geometrical information about the throat 
in the compact internal space. Current proposals for the 



form of /(</>) include f(<p) 



(a constant), for the 



case of an AdS throat. Another case considered in the 
literature is a constant function /(</>) = a [8l|. In this pa- 
per, we will keep f{4>) and V{<p) as general non-negative 
functions, i.e. V(<fi) > and /(</>) > 0. 

In a spatially flat Friedmann-Robertson- Walker 
(FRW) metric with scale factor a(t), X = <j) 2 /2 for a 
homogeneous field 4>, and it can be shown that the en- 
ergy density and pressure of the DBI field are given by 



7 + 1 



4> 2 + V(4>), 



P4> = ^—^-V{4>), 

7 + 1 



(3) 



where a dot denotes a derivative with respect to the cos- 
mic time, t, and 



f 



i - f(tW 



(4) 



which characterises the motion of the brane and serves 
as the Lorentz factor. As in usual special relativity 7 > 1 
because /(</>) is non-negative. 

If we also take into account a perfect fluid with equa- 
tion of state p m = w m p m (—1 < w m < I), the basic 
cosmological background equations are given by 



H 2 = l 
3 

•• 3H 

</>+—< 

T 



7 + 1 



4> 2 + V(<(>) + Pn 



Vj, /,* (7 + 2)(7- l)ig 
7 3 2/ (7 + 1)7 



2 = O, 



3if(l + w m )p m = 0, 



(5) 



(6) 
(7) 



where the subscript <f> means differentiation with respect 
to the field, H = - is the Hubble parameter and we have 
set 8nG = 1 for simplicity. Eqs. (|5|) - (J7J close the system 
that determines the dynamics. 

In order to obtain the late-time attractor behaviour of 
the system we introduce the following set of dynamical 
variables: 



7^ 



V 



V 



V3H' 



7 



1 

1 



(8) 



By defining x and y in this way it is clear that in the 
limit of 7 — > 1, we recover the dynamical variables for a 
canonical scalar field originally proposed in Ref. [57[ . By 
doing this, we can intuitively think of x corresponding 
to the contribution of the kinetic energy, and y to the 
potential energy of the field. For another dynamical de- 
gree of freedom which arises due to the introduction of 
the function f(<fi) in this set-up, we adopt 7 as our third 
dynamical variable. This particular choice is intended to 
make the phase space being compact. 

In terms of these variables, the Friedmann constraint 
© can be expressed as 



x 2 + y 2 + Q m = 1, 



(9) 



where Sl m = p m /(3H 2 ). The energy fraction and equa- 
tion of state of the DBI field are given by 



= x 2 + y 2 , 



~ 1 1 
jx — y 

x 2 + y 2 



(10) 
(11) 



Similarly, for a general potential V(<j>) and brane ten- 
tion /(</>), by introducing A and p which are defined by 



A 



V 



(12) 
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Eqs. ([6]) and ([7]) can be written as 



A. Autonomous System 



-v/37(l+7)Ay 2 

+~x[(l + j)(x 2 - 1) + (1 + w m )(l - x 2 - y% 

(13) 



y,N = 



~f,N 



+ 7> 2 + (1 + - z 2 - 2/ 2 )], (14) 

7(1 -I 2 ) 



3^1 + 7 + V3^-(M^ 2 - Ay 2 ) 



V37(l+7)A 2 (r-l)z, 
V37(l+7)M 2 (S-l)x : 



(15) 
(16) 



In order to represent Eqs. ©-([7} by an autonomous 
set of equations, we first introduce the variables A and jl 
denned by 



A= - 



V, 



fqVi +1 



fq+^yq ' 



(19) 



where q = —l/(p + r). Initially we assume p + r =/= 
0, and will address that special case later. Given the 
definitions in Eqn. (|19[) . we can easily verify that A and 
jl are constants given by 



A 



aivi' 



-er 



-A, 
V 



(20) 



where N = In a, x,n 
r and S through 



dx 
tlN 



etc... and we have introduced 



VV, 



//,, 



V 



2 ' 



(17) 



Although Eqs. (|13|) - (|T6)) hold for a general poten- 
tial and brane tension, if both are given as exponential 
functions of <j), V oc exp[— \<fi], f oc exp[— fjuf>], then from 
Eqs. (fT5|) and (fT6]). A and (i become constant and Eqs. ([T3|) 
- (fT5j) form the closed autonomous system. We will see 
the late-time attractor structure of this case in Sec. IIVI 

There is an important new feature which emerges in 
the DBI case and is not present for the case of the canon- 
ical scalar field. In the latter case only an exponential 
potential can truly lead to an autonomous system, but 
as we shortly show, this fact does not apply to the DBI 
field. By combining the degree of freedoms related with 
V(<j)) and f((p), we will be able to obtain an autonomous 
system for a wider class of potentials and brane tensions 
(see also (69|). 



III. MODELS WITH A GENERAL POWER-LAW 
POTENTIAL AND BRANE TENSION 

Here, we consider models where the scalar potential 
and the brane tension are arbitrary non-negative power- 
law functions of the DBI field 



V{fl>)=aW, f{<l>) = v\W, 



(18) 



with constants a > and v > 0. 

Guo and Ohta ,69] analysed such a model for the case 
p = 2 and r = —4 but as we will now show it can be ad- 
dressed as an autonomous system for more general com- 
binations (the exception being r + p = which requires 
a separate treatment as we will see). 



where e is 1 for cj> > and is —1 when < 0. The 
requirement that A > is equivalent to demanding < 
for p > 0, and <p > for p < 0. We, therefore, restrict 
our analysis to A > without loss of generality. Then, 
for physically interesting cases, we can also restrict our 
solutions to those where <f> > (x > 0). 

From the form of V(^>) and /(</>) in Eq. (fT%|) it follows 
that A and jl are related with A and fi through 



A = 



7(1 + j)x 2 



M = 



(l-7)y 2 
7(1 + 7)x 2 



(21) 



Therefore since X.n and fi t w can then be expressed in 
terms of x,n, d.n, 1,n, defined in Eqs. (fl~3"|) - (|16[) . we 
only need to solve those three equations. 

In terms of the dynamical variables we have defined 
above, the basic equations can be expressed as 



X.N 



V3(l + 7 )t 



A7*-«(l - 7)9 



29+2 



r 2q 



(l+ 7 )(a; 2 -l) + (l+i i ; m )(l-a; 2 



y,N 



(22) 



+ -y[(l + 7 )x 2 + (1 + w m )(l -x 2 ~ y 2 )] , (23) 



1,n = 7(1 - 7 ) X 



3- , ^A75" 9 (l-7) 9 r ' 1 ' 



1 + 7 



r 2 9 +l 



2 2 

-x - y 



(24) 



Equipped with the basic equations forming an au- 
tonomous system, we will peform the stability analysis 
to obtain the late-time attractor structure in the follow- 
ing. 
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B. Standard Fixed-point Solutions 

Here we obtain fixed-point solutions of the dynami- 
cal system given by Eqs. (j2"2")) - ([2"4")) . The fact that there 
can be terms involving x, y and 7 in the denominators 
of these equations means that we have to treat the cases 
where these terms vanish with considerable care. In this 
subsection we make sure we are working in regimes where 
any possible ambiguities involving possible divisions by 
zero or ratios of '0/0' are avoided. In order to distinguish 
between them, we call the fixed-point solutions consid- 
ered in this subsection to be 'standard' fixed-point so- 
lutions. We will address the other cases in subsection 

From Eq. (|2~4"j) . the requirement that 7 is a constant at 
the standard fixed points implies three possible scenarios: 

(a) 7 = 0, 

(b) 1 = 1, 

w s =v^^" f(1 -^(-; aS+ ^ 

(25) 

We will investigate each scenario in turn in what fol- 
lows. 



1. Case a : 7 = 

This result is only valid as long as q < 1/2. In this case, 
since the Lorentz factor ((4]) tends to infinity, we shall 
refer to this class of fixed-points as the "ultra-relativistic" 
solutions. Substituting 7 = into Eqs. ([2"2" 1) -(|2"3 ]) and 
eliminating the terms that obviously vanish we obtain 



V3r-i y 2q+2 

+ ^x[(x 2 - 1) + (1 + w m )(l - x 2 - y 2 )] ,(26) 

4- 1 a y 2 i +1 
= --^^ 

+ 3 -y[x 2 + (l + w m ){l-x 2 -y 2 )], (27) 



1,N 



2<; 



2<j+l 



P 



y 



(28) 



Due to the freedom in choosing the sign and the mag- 
nitude of q, it is clear that there are situations where 
some of the terms in these equations are ill-defined. We, 
therefore, classify our analysis in terms of the value of q. 

For < q < 1/2, we find the following standard fixed 
points: 



(01) (as,y,7) = (1,0,0) 

(02) (x,2/, 7 ) = (*„, 0,0) 

with < xq < 1, w„ 



(29) 
(30) 
(31) 



We refer to the point (al) as the ultra-relativistic ki- 
netic dominated solution, and to the family of solutions 
(one for each :ro) (a2) as the ultra-relativistic kinetic-fluid 
scaling solution. For q < —1/2, we find the following 
standard fixed-point solution 



(o3) (x,y,j) = (0,1,0) 



(32) 



which represents the ultra-relativistic potential domi- 
nated solutions. 

There are no 7 = solutions for the case, —1/2 < q < 
0, but for the special case of q = 1/2, there exist two 
other standard fixed points in the system of equations 



(o4) (x,y,j) 



X(V\ 2 + 12 — A) -A + Vl2 + A 2 



2^3 



(o5) (x,y,i) 



3(1 + w m ) 3 \/3(l + w m ) 
\ 2 w m A 

y/3(l + w m ) 



with w m < 0, A > 



\J-w r , 



(33) 



which we call the ultra-relativistic kinetic-potential scal- 
ing solutions, and the ultra-relativistic kinetic-potential- 
fluid scaling solutions, respectively. 

It is worth commenting on the solutions (al) — (ah). 
(al) is well known and corresponds to the case where the 
kinetic energy of the DBI field dominates over the po- 
tential energy, leading to an effective equation of state 
for the DBI field in Eq. (ITT]) which mimics that of dust 
(iVcf, = 0). The existence of the solution (a2) where this 
ultrarelativistic kinetic dominated DBI field scales with 
matter was previously obtained by Ahn. et. al [13, H3 
who also obtained the solutions (a3) and (a4). In fact 
(a3) is the solution actively investigated in the context 
of DBI inflation, for example as seen in Refs. The 
solution (a5) is of a new type, where the ratio of the 
kinetic and potential terms remain constant in the ultra- 
relativistic limit. The particular case of this solution for 
p = 2 and r = —4 was first discovered in 69], although 
here we have shown that this type of solution exists as 
long as p + r = —2. 



2. Case 6:7 = 1 

In this case the Lorentz factor 7 = 1 and the DBI 
field mimics the behaviour of a canonical scalar field. We 
shall refer to this class of fixed-points as the "standard" 
solutions. Substituting 7 = 1 into Eqs. (|2"2"j) - (j2~4"l) . we 
obtain 



5 



= — A(1 ~ 7) 

+ ^x[2(x 2 - 1) + (1 + w m )(l -x 2 ~ y 2 )] 



(34) 



VQ~ v 2q+1 



+-y[2x 2 + (I + w m )(l - x 2 - y 2 )] , (35) 



7* = VeMi-ir^^-v 2 



(36) 



Following our previous approach we classify the stan- 
dard fixed points of this system based on the range of 
values the parameter q can take. We find that for q > 0, 
there exists a standard fixed point 



(61) (x,y,j) = (1,0,1) , 



(37) 



which is the standard kinetic energy dominated solution. 

For the special case of q = (p + r — > oo), we find two 
extra interesting standard fixed-points: 



<> 2 > f^^ 1 , 

with A < V6, 



(38) 



and 



with A > \/3(l + io m ) and — 1 < iu m < 1, 

(39) 

which we call the standard kinetic-potential scaling so- 
lutions, and the standard kinetic-potential-fluid scaling 
solutions, respectively. 

In reviewing these solutions, recall that the DBI field 
behaves just as the usual canonical scalar field when 
7 = 1, hence the fixed point solutions found in Case 
(6) are already well known. For example the properties 
of (62) and (63) are identical with that of a standard 
power-law inflationary solution and scaling solution re- 
spectively, obtained with an exponential potential in the 
presence of a canonical scalar field |53l - l57| . In fact (62) 
was obtained previously in the context of DBI Inflation 
in M M. 



Case c : 3 



Here 7 is a constant which is different from or 1. In 
this case, since the Lorentz factor 7 > 1 and is constant 



as defined in Q , we shall refer to this class of fixed point 
solutions as the "relativistic" ones. 

As there are no values of q which permit standard 
fixed-point solutions x = y = in equations (|2"2"|) - (|2"4"1) . 
we therefore begin by exploring the possibilities of either 
x or y being in this case. 

For the case of x ^ 0, y — 0, we find the following 
standard fixed points in the system which exists only for 
q = 0: 



(cl) (x,»,7)= 1,0,- 



and 



(c2) (x, y, 7) 




(40) 



With < W m < 1, ,i < _ VW+Wm) ^ (41) 



where we refer to (cl) as the relativistic kinetic domi- 
nated solution, and to (c2) as the relativistic kinetic-fluid 
scaling solution. 

For the case x — 0, y ^ 0, a standard fixed-point 
solution exists but only for q = —1/2 



(c3) (x,y,j) 



0,1, 



V3 



Vx 2 



(42) 



which we call the relativistic potential dominated late- 
time solution. 

For the remaining cases with x 7^ 0, y 7^ it follows 
that combining Eqs. (|2"2"|) - (|2"31 for general q, yields the 
condition 



3 = 



3 -A 7 3-9(l-7)? 



1+7 



Comparing this with Eq. (f2~5j) , we see that in Case (c) 
the condition for a fixed point for general non-zero x and 
y requires p = — r (or q — > —00) a limit we have decided 
not analyse in this section. 

Summarising the results of Case (c) we note that the 
fixed-point solutions (cl) and (c2) are completely new, 
while (c3) can also be derived as a special case of line 
5 of Table I in (52[. Of particular note for cosmology 
is (c3) which is the concrete example of an inflationary 
solution with constant 7 which differs from 1 and 0. Phe- 
nomenologically this is a very interesting solution when 
considering cosmological perturbations in DBI inflation. 

In TABLE HI we have provided a breakdown of the 11 
standard fixed-point solutions obtained in these class of 
models corresponding to cases (a) - (c). 
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C. Stability Analysis for standard fixed-points 

We now turn our attention to carrying out a stabil- 
ity analysis for the standard fixed-points obtained in the 
previous section. Calling these points in general x c , y c 
and 7c we consider small fluctuations about them given 
by 



Sx , 



V = Vc + Sy . 



7 = % + Sj, 



(44) 



and consider solutions of the form 8x oc e^^, Sy oc e wN 
and <$7 oc e wN . As before we consider each case in turn 
and examine the dynamical behaviour of the system close 
to their fixed-point positions on the phase plane. 



1. Case a : 7 = 

Expanding Eqs. (|22|) - ([24"1) around the ultra- relativistic 
kinetic energy dominated solution (al), the ultra- 
relativistic kinetic-fluid scaling solutions (a2) and the 
ultra-relativistic potential dominated solution (a3) we 
obtain the following eigenvalues 



(al) : wi = 3 , w 2 = - , w 3 = -3u>„ 
3 

(o2) : wx = 3 , w 2 = - , w 3 = 



(45) 
(46) 



(a3) 



2' w 2 =-3(l + w m ), w 3 = 3.(47) 



The presence of positive eigenvalues in each case indicates 
that all three solutions are always unstable. 

For the ultra-relativistic kinetic-potential scaling solu- 
tion (<x4), we obtain the following eigenvalues 



W\,2 



3(4w m + 8 + A 2 - Vl2 + A 2 ) 



1 



w 3 



-\+yi2 + X 2 j (A + /i) 

where A(X, w m ) = 8 ^-12 - A 2 + A \/ 12 + A 2 
-6w m — 6 — A 2 + Wl2 + A 2 



9(-A 2 + AV12 + A 2 - 4(2 + w m )) 2 



(48) 



It can be seen that wx t 2 are non-positive if the con- 
dition w m > — 1 + (A/6)(\/A 2 + 12 — A) is satisfied. 
Rewriting we see that this condition becomes < A < 
V3(l + w m )/y/ —w m if w m < 0. Therefore, considering 



also W3 we see that these solutions are stable if the addi- 
tional condition X + p, < is also satisfied. However from 
Eq. (f2T)|) we know that A + fi = ^ (r + p) and given that 
the condition for the existence of the solution to (a4) is 
q = 1/2, we see that this solution is stable only for the 
case p > where r + p = —2. 

A similar analysis of the ultra-relativistic kinetic- 
potential-fluid scaling solution (a5) yields 



Wx,2 = ^(1 - Wm) 



l±Jl + B(X,w m ) 



w 3 = 3(1 + w m ) 



X + p 



V, K(\ \ 8(^mA 2 +3(l + w m ) 2 ) 

where B(X,w m ) = . (49) 

X 2 w m (l + w m ) 

In this case w m < and B < for A > "\/3(l + 
w m )/y/—w m , hence wx,2 are negative whenever these so- 
lutions exist. Of course the bound on A is the comple- 
ment of that arising above for the stability of (o4) , hence 
it follows that only one of the two solutions (a4) or (a5) 
can be stable for a given value of A. The other condition 
required for stability is once again that w 3 < 0, which as 
before corresponds to, A + p < or p > (recall we are 
assuming A > 0). It is worth noting that A = \p\j ^fav 
for q= 1/2. 



2. Case 6:7 = 1 

For the standard kinetic energy dominated solution 
(61) the eigenvalues are 

wx = 3(1 — w m ) , w 2 =3, W3 = — 6, (50) 

for q ^ and 

wx = 3(1 - w m ) , W2 = 3-^A, w 3 = -6 — V&p, 

(51) 

for o = 0, which clearly indicates that these are unsta- 
ble solutions. For the standard kinetic-potential scaling 
solutions (62), (also q — 0), we find the eigenvalues 

A 2 - 6 

wx= — - — , w 2 = A 2 - 3(1 + w m ), w 3 = ~X{X + p), 

(52) 

s uggesting t hese solutions are stable when A < 
\/3(l + w m ) and A + p > 0. It is worth noting that 
A = \p\ and A + p = — \p\/(pq) for q = 0. Then, the 
condition A + p > corresponds to p < for q — > +0 and 
p > for q ->• -0. 

Similarly, for the standard kinetic-potential-fluid scal- 
ing solutions (63) we obtain 



Wl,2 = |(1 - Wm) 



7 



-1±4/1 



8(3(1 + w m ) - A 2 )(w m - 1)(1 + w m ) 



\ 2 (-l+w m f 



w 3 = -3(1 + w m ) 



A + /2 
A 



(53) 



Since A > y3(l + u> m ) whenever these solutions exist, 
and W2 are always negative. Therefore, these solutions 
are stable if the condition X + fl > is satisfied. Note also 
from the condition above for the stability of (62), if the 
conditions for both (62) and (63) occur, only (63) will be 
the stable late time attractor, in other words the standard 
kinctic-potential-fluid scaling solution will dominate over 
the standard kinetic-potential scaling solution. 



3. Case c : 3 



7) 9 ^r (- 



;x 2 + y 2 ) 



it is dominated by the kinetic term. It is also worth 
noting that since 7 = 3/(/2 2 — 3) for this fixed-point, the 
stability conditions are given in terms of p, as A 2 > fi 2 
and w m > 3/(/i 2 — 3). We can go a little further using 
Eq. (|20|) implying A = — ep. It follows that the condition 
for stability is 

have the opposite sign whilst satisfying q = 

For the relativistic kinetic-fluid scaling solutions (c2), 
a similar analysis produces 



«?i,2 = |(1 - w m ) 



< 1 with ^ < 0, so r and p have to 



-1± a/1 



8(3(1 + w m )(p, 2 w m - 3(1 + w m ))) 



+w m ) 



w 3 



3(A + +w m ) 
2/2 



(55) 



For the relativistic kinetic energy dominated solution 
(cl), expanding around this fixed point yields the folow- 
ing eigenvalues 

wi = 3(-w m + 7) 



w 2 = 



w 3 



3 + 37-\/3A^7(l + 7 ) 
2 

3(-l + 7) 



(54) 



where W3 is clearly always negative. The parameter space 
for which this solution is stable is A > \/3(l + 7V7 and 
w m > 7- This is natural since if the potential is steep 
enough, the kinetic term easily dominates the potential 
term and if w rn > j(= w^), the energy density of the fluid 
decreases faster than that of the DBI field even though 



Since p? > 3(1 + w m )/w m for these solutions, either w\ 
or u>2 is always positive, which means these solutions are 
unstable. 

The relativistic potential dominated solution (c3) 
yields the following eigenvalues 

twi=0, w 2 = -3, w 3 = -3(1 + w m ) . (56) 
As this has a zero eigenvalue, we say this is a 'marginally 
stable' solution in the sense that there is no instability 
growing exponentially, although it could be unstable to 
higher orders in the perturbation. Obviously, the stabil- 
ity of this point is weaker than that of the fixed-point 
with three negative eigenvalues. 

To help clarify all the possible standard fixed-point 
solutions and their stability the information just provided 
is summarised in Table U below. 



D. Fixed-points arising from solutions that initially 
appear singular 

As mentioned earlier, for models with a general power- 
law potential and warp factor, in addition to the usual 
standard fixed points, it is necessary to check if the points 
where the denominator is in the right hand side of 
Eqs. (|2"2"1) - ([2~4")1 can also be late-time attractor solutions 
for the system, a situation which does not arise in the 
case with a canonical scalar field 57]. 

Even in the case that the point itself leads to a singular- 
ity, it doesn't necessarily mean the system is ill defined. 
For example it could be that the solutions approach the 
point exponentially slowly (i.e. like exp[— N)), hence it 
would take an infinite time to reach the singularity and 
physically there is no ill behaviour in the system. In par- 
ticular as long as the ratio of the singular terms (loosely 



called '0/0') tends to a constant value then the system 
can be analyzed for the stability of these fixed points. As 
standard techniques can be applied to judge the stability 
of such a point, we also call them fixed-points in what 
follows. In the following, depending on the value of q we 
show there are 6 kinds of fixed-point where '0/0' is finite 
in the phase space. As in the previous analysis, we have 
excluded the special case with p = —r here. 



1. Standard potential dominated solutions 

First, we consider the point (al):(a;, y, 7) = (0,1,1) 
which in the case of a canonical field would simply cor- 
responds to a standard slow-roll inflationary solution. 
However, in this case x = leads to a singularity in 
Eqs. (J22J)-(E3J) for q > -1/2 and 7 = 1 does the same for 
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a; 


J/ 


7 




Valid g 


Existance 


Stability 


(al) 


1 








1 


< q < 1/2 


VlU m ,VA, V/i 


unstable 


(o2) 


zo 








Kg 


< <7 < 1/2 


< s < 1, ui m = 0, VA,V/i 


unstable 


(o3) 





1 





1 


<? < -1/2 


Vw m ,VA, V/i 


unstable 


(o4) 


/ A(\/a 2 + 12 — A) 







1 


g = 1/2 


Viu m , VA, V/i 


Wm > — 1 + x\, P > 


V 6 




(o5) 


/ 3(l + u,,„) 3 
V ^ 2 ">- 


VS(l+w m ) 
A 





3(l+iu m ) 2 
A 2 ™ m 


g = l/2 


Wm <0,A> ^ 1 +""" ) ,va 


p ^ y3(l+m m ) 

to m < 0, p > 


(61) 


1 





1 


1 


q > 


VTO m ,VA, V/i 


unstable 


(62) 


A 

76 


v D — " 

V6 


1 


1 


g = 


V™ m , A < \A5, V/i 


\P\ < V 3 (1 + W m), 

l + ^>0 


(63) 


\/2\ 


V3(l-^) 
V^A 


1 


A 2 


g = 


Vw m ,A> ^/3(1 + «;,„), V/i 


P > \/ 3 (l + Wm )> 

! + i >° 


(cl) 


1 





3 

M 2 -3 


1 


g = 


/i < — \/6, Vu> m ,VA 


> ^-tt(> 0), 

o>^>-i 


( C 2) 


^/3(l+w m ) 





w m 


3(l+t» m ) 


g = 


0<«, m <l,A< v ' 3C ^ m) ,VA 


unstable 


(c3) 





1 




1 


9 = -1/2 


VTO m ,VA, V/i 


marginally stable 





TABLE I: Summary of the standard fixed-points and their stability in the models where the potential and brane tension are 
arbitrary non-negative power-law functions of the DBI field. Notice that we have restricted A > and x > 0. We have expressed 
the conditions for the stability in terms of (p, r) not (A, /i) as these are the more fundamental quantities. 



q < 0, hence we have to tread carefully in analysing the 
system. 

Since the coordinate y = 1 does not result in any ill- 
defined behaviour in Eqs. (f22|) -(|24 |) . we can consider the 
reduced system in which we determine the leading order 
behaviour of x and 7 around the point (al). Writing 
x = Q + 5x, j = I — Sj and keeping only the leading 
order terms, we obtain 



Sx M = 



N 



(57) 
(58) 



If we introduce j3 



-, for q < -1/2 and q > 0, 



this appears to be of the form of '0/0' and requires a 
careful analysis to properly understand the behaviour of 
this system. In this case from Eqs. (|57p - (1551) . it becomes 



(Sx) 



2q+l 



-3N 



(J7)« 



41. 

6 



(59) 



where c is an integration constant. At late times (N — > 
00) we see that even though both terms tend to zero, 
the ratio (Sj) q / (8x) 2q+1 approaches the contant given 
by (Sj) q I (Sx) 2q+1 = \/6/A. Then, it is possible to sensi- 
bly discuss whether the point (al) = (0, 1, 1) is stable or 
not in terms of the remaining two-dimensional system ob- 
tained by substituting (1 — j) 9 /x 2g+1 — v6/A back into 
Eqs. ([2"2"j) - (j2"3")l . It is worth noting that as the coefficient 
of N in the exponential function in Eq. (I59[) is —3 which 



is negative, this solution is stable along the direction of 

The eigenvalues w±, W2 and W3 corresponding to evo- 
lution of the perturbations in x, y and 7 respectively are 
obtained from Eqs. ([57 J) .(|2"3" ]l and JST 



(60) 



wi = = w 3 , 



w 2 = -3(1 + w m ) . 



To be specific the results w\ = = W3 arise be- 
cause at leading order, the analytic solution including 
(8"f) q / (5x) 2q+1 — \/6/A leads to a vanishing right hand 
side for Eqs. (1571) and (|55|) . The fact that at leading 
order all eigenvalues are non-positive, implies that for 
q < —1/2 and q > there are solutions which tend to 
(al) regardless of the values of A, jl and w m even though 
it is not a fixed-point solution in a usual sense. 

But this zero eigenvalue implies that the stability of 
this point is weaker than the point which has negative 
eigenvalues for all three directions. As in the case (c3), 
we say this solution is 'marginally stable'. Of course, to 
obtain the strict stability of this solution we would have 
to go to higher order. 



2. Ultra-relativistic potential dominated solutions 

Next, we consider the point (a2):(x, y, 7) = (0,1,0) 
whose behaviour is the same as (a3) except for the fact 
that x = is singular for q > —1/2 and 7 = for q > 1/2 
in Eqs. ([22])-([2l]>. 

Following the arguments used for (al) and writing x = 
+ Sx, 7 = + #7 we obtain to leading order: 



9 



Sx 



N 



Sj,N 




(61) 
(62) 



For -1/2 < q < 1/2, p = 



once again appears 
ill denned, however we can solve Eqs. (1BT1) - (j6"2")l to give 



(Sx) 



2q+l 



(50 



ce 



-3N 



3 



(63) 



where c is an integration constant. At late times (N —> 
oo ) we see that even though both terms tend to zero, the 
ratio (Sj) 1 ^-i/(Sx) 2 i+ 1 -> V3/A. 

As in the case of (al), it is possible to discuss whether 
the point (a.2) = (0,1,0) is stable or not in terms of 
the remaining two-dimensional system obtained by sub- 
stituting the particular solution ^ 1 / 2 ~i / x 2q+1 = \/3/\ 
back into Eqs. ([22"]) -(j23" |) . It is worth noting that as the 
coefficient of N in the exponential function in Eq. (|6"3")l 
is —3 (i.e. negative), this solution is stable along the 
direction of (3. 

The corresponding eigenvalues for the perturbations in 
x,y and 7 follow from Eqs. (pjlj) . (|2"3"|) and (p32"j) and are 
given by 



wi = = W3 . 



w 2 = -3(1 + w m ) . 



(64) 



Note that the eigenvalues for (a2) are identical to those 
of (al) hence to obtain the full stability of the system, 
as in that case, we would have to go to higher order to 
obtain the strict stability of the system 



3. Ultra-relativistic kinetic dominated solutions 

Next, we consider the point (aS):(x, y,j) = (1,0,0) 
whose property is the same as (al) except that y = is 
singular for q < and 7 = is singular for q > 1/2 in 
Eqs. ©-d!. 

In this case, since the coordinate x = 1 does not re- 
sult in any ill-defined behaviour in Eqs. (|2"2"]) - ([2"4|) . we can 
consider the reduced system in which we determine the 
leading order behaviour of y and 7 around the point (a3). 
Writing y — + Sy, 7 = + £7 and keeping only leading 
order terms, we obtain 



Sy,N 

Sj,N 



(^, 



V p ( s v) 9 



(65) 
(66) 



Following the earlier examples, for q < and q > 1/2, 



we introduce (3 



(«57)V2 



and proceed to show that 



it has a well behaved non-trivial behaviour of its own. 
Solving Eqs. ([6"5j) -([6l) |) we find 

(5y)- 2q _ 



-#7V 



(50 



\/3~ 



2q + (2q 



D- 

p 



(67) 



where c is an integration constant. At late times (N ~ > 
00) we see that even both terms tend to zero, the ra- 
tio (5~f) 1 / 2 ~ q / (5y)~ 2q approaches the constant given by 
^i/2-9/(fy)-29 = v /3/[A(2 ? +(2g-l)r/p)]. Since both 
6y and ^7 are positive, for this solution to be physical, 
for a given q, r/p should satisfy 



2q + (2q - 



r 

l)->0, 

p 



(68) 

-2 or p < with 



which is equivalently p > with r < 
r > -2. 

As in the previous cases, it is possible to discuss 
whether the point (a3) = (1,0,0) is stable or not in 
terms of the remaining two-dimensional system obtained 
by substituting 7 1 /2-?/ y -2g = ^3/[A(2g + (2g - l)r/p)} 
back into Eqs. (|22l) - (|2"3"|) . It is worth noting that as the 
coefficient of N in the exponential function in Eq. (|67|) is 
—3/2 which is negative, this solution is stable along the 
direction of j3. 

The eigenvalues for the perturbations in x, y and 7 
follow from Eqs. fl55J} , flBU) and (J66j) to give 



Wl 



W2 = 



—3w„ 
3(2q 



p' 



2[2q- 



(29-1)5] 



6g(l 



w 3 



[2g + (2g-l)j 



2 + r' 



(69) 
(70) 

(71) 



where in Eqs. (I70|) - (l7l1) we have made use of the fact 
q = —-pj^- Although W2 and W3 diverge for r = —2, 
from Eq. (|B"5)) this case is not physical. 

From these eigenvalues, we see that for both cases q < 
and q > 1/2, stable solutions require w m > and 
r < —2. Note that in Eq. (|68|) this implies p > for 
stable solutions. For q < 0, the additional constraints are 
r/p > —1, whereas for the case q > 1/2 it is — (2/p+l) < 
r/p < —1. 



4- Ultra-relativistic kinetic-fluid scaling solutions 

In this case the behaviour is identical to (a3), for the 
special case with w m = 0. However, there is another 
important point (a4):(x, y, 7) = (xq,0,0) with Xq satis- 
fying < xq < 1 and whose property is the same as (a2) 
in that there are a family of solutions depending on the 
value of xq chosen. Again y = appears singular for 
q < and 7 = for q > 1/2 in Eqs. p^ - (|23|) . 

Following the arguments used for (a3) and writing x = 
x o, U = + Sy, 7 = + S'f, we obtain to leading order 



Sy.N 



V3 



Sy \/3 - \x 



l-2q 



(Sy) 



-2q 



(72) 



10 



fry n = V3S"f 



\/3 + -Axi 



-2q 



my- 



(Sy) 



-2q 



(73) 



For q < and q > 1/2, again introducing ft 



my 



we can solve Eqns. (|72|) - (|73|) to give 



(Sy) 



-2q 



(57)1/2-9 



2« 



2< 7 +(2< 7 -l)- 



(74) 

where c is an integration constant. Once again at late 
times (N — > oo) we see that even though both terms 
tend to zero, the ratio (d)^) 1 / 2-9 / '(Sy)~ 2q approaches the 
constant given by (Sj)^ 2 - q /(5y)- 2q = V?>/[\x 1 - 2q {2q + 
(2q — l)r/p)]. Since both Sy and #7 are positive, for this 
solution to be physical, for a given q, r/p should satisfy 
Eq. ([6"5)> . It is worth noting that the coefficient of N in 
the exponential function in Eq. (1741 is —3/2 which being 
negative implies this solution is stable along the direction 
of ft. 

The eigenvalues for the perturbations in x, y and 7 
turn out to be identical to those for (a3) in Eqs. f7D|) 
and (171]) . except for the case of w\. In particular from 
Eqs. ([22J) , (J72l) and J73|) we obtain 



(75) 
(76) 



w 3 



3(2g -!)(! + ^) _ 3 
2[2q+(2q-iy-} 2 



[2« + (2g-l)£] 2 + r' 



(77) 



Again, although W2 and W3 diverge for r = — 2, this case 
is excluded as ft is then no longer a finite constant. 

wi shows that although there is a zero eigenvalue along 
the x direction, when w 2 < 0, there are solutions ap- 
proaching the point (a4). From the discussions for (a3), 
we see that for both cases q < and q > 1/2, stable so- 
lutions require r < —2, p > 0. For q < 0, the additional 
constraints are r/p > —1, whereas for the case q > 1/2 
it is -(2/p + 1) < r/p < -1. 



5. Standard fluid dominated solutions 

Next, we consider the point (a5):(x, y, 7) = (0,0,1) 
which in the case of a canonical field corresponds to the 
usual fluid dominated solution. However, here x — is 
singular for q > —1/2, y — and 7 = 1 are singular 
for q < in Eqs. ([2"2 ]) -(|2"l )) , Since all the coordinates 
x = 0,y = 0,7=l can result in ill defined behaviour 
in Eqs. (f2"2"|) - (|24[) . the stability analyis for this point is 
more complicated than in the previous examples in this 
subsection. 

To investigate the behaviour of perturbations in x, y 
and 7 around (a5), we write x = + Sx, y = + Sy, 
7 = 1 — ^7 and keep only the leading order terms, to 
obtain 



Sx 



Sy. 



N 



N 



Sx 
~2 

'111 
2 



. 3(1 _„ m) + V 5A( fa '" I - 2 '»' W 



3(l + m m )-V6A 



(*7)-« 

(fa)- 2 ^ 1 ' 



(Sy)- 2q (S^)- q 



(78) 
(79) 



S~f,N = —257 



2"((5?/)- 2 «(57)-9 Vp 



(80) 



We can in fact make considerable progress by solving 
for combinations of the variables as in the previous ex- 
ample. Introducing 6 = ^ ) x X q ^-y q 
Eqns. ([78 ]) -([80 l we obtain 



and P = ^3-, from 



N 



-(1 - m m ) + 



-2q 1 



^ = 6-^6^(1 + ^), 



(81) 
(82) 



These equations pick up two of the three degrees of 
freedom of the system given by Eqs. (|78)l -(|80|). and sta- 
bility in this two-dimensional system is necessary for the 
stability in the full three-dimensional system. We find 
this two-dimensional system has three fixed-points char- 
acterised by 



?,/3) = (0,0) , and 



3(1 - w m ) 
-2qV6\ (l + 



and 



(83) 
(84) 



3(1 



j 4q(l + ^+(l-w m y 



V6A 



2q 1 



(85) 

It turns out that the eigenvalues corresponding to 
fluctuations about Eq. (153")) are —3(1 — w m )/2 and 6 
whilst those corresponding to Eq. ([51)) are 3(1 — w m ) 
and 3[(1 — w m ) + 4q(l + r/p)]/[2q(l + r/p)}, indicating 
that both solutions are unstable. 

Turning to the stability of the fixed-point correspond- 
ing to Eq. ([55)) we first of all note that in order for the 
solutions to be physical requires a number of conditions 
be satisfied. For 9 and ft to be finite constants requires 
q < since Sx and Sy have the same time dependence. 
Further, since Sx, Sy, S^ are all positive, 9 and ft must 
be positive, which leads to the conditions 



(86) 



iq[ 1 + - ] +(l-«j m ) >0, (87) 



1 + 2q [ 1 + - ) > . 
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which are equivalently p > 4/(1 — w m ) or p < 0. 

Now, by considering small perturbations around this 
point of the form 59 = Ae w « N + Be w " N , 8(5 = Ce w - N + 
De WbN , we obtain the following eigenvalues, 



w a ,b 



3(2q{3 + w m )(l + l) + (l-w m )) 



4(l + 2g(l + £)) 



-1 ± 



1 - C(A, W m , -) 
P . 



where C(X, w m , —) is 
P 

8(1 + 2q(l + + Wm )(4g(l + £) + (1 - Wro )) 

(2<?(3 + Wm )(l + £) + (l- Wm )) 2 



and (+,-) is associated with (iy a , itfb) respectively. Given 
the conditions (|86|) and ([87]) are satisfied for physi- 
cally relevant situations, it follows that stable eigenvalues 
(wa,b < ) are obtained if 



2q(3 + w m )(l + 



P 



(1 



> 0. 



(89) 



which is equivalent to p > 2(3 + w m )/(l — w rn ) or p < 0. 

In fact ([5"9"|) is more restrictive than ([55]) and ([57]) as 
long as we are in the physically acceptable range — 1 < 
w m < 1. 

If the condition (|89l) is satisfied, the terms including Sx, 
Sy, Sj in Eqs. (1751 - ([50]) are all finite and constant. Then 
it is possible to sensibly discuss whether the point (a.5) = 
(0, 0, 1) is stable or not in terms of the remaining one- 
dimensional system obtained by substituding ([55]) back 
into say Eq. ([2"2]). 

Then we find the eigenvalue for the perturbations in x 

3q(l + w m )(l+ r -) 3(1 + Wm) 



l + 2 g (l + J) 



P 



Notice that we have already obtained the stability for 9 
and {3 directions, the stability based on the eigenvalues 
for the directions of x, y and 7 have the same infor- 
mation. Actually, the eigenvalues corresponding to the 
perturbation in y and 7 are given by 



U'2 



W 3 = 



3g(l + Wm )(l + p 

l + 2g(l + §) 
3(l + Wro )(l + £) 
l + 2g(l + §) 



3(1 + w m ) 



(91) 
(92) 



which are consistent with 



wi = w 2 , (93) 
(-2q + 1)W! + 2qw 2 + qw 3 =0, (94) 

derived from 9 and (3 are constant. Although w\, w 2 and 
u>3 diverge for p = 2, from Eq. ([59"]) . this case is excluded. 



We can summarise the constraints on the allowed pa- 
rameters if the solution (a5) is to be stable. From 
q < we have r + p > and from w± < we have 
p > 2 and 1 + r/p > 0. Then, from Eq. ([55]). we have 
p > 2(3 + UJ m )/(l - iu m ) > 2 and w m < (p - 6)/(p + 2). 



ff. Ultrarelativistic fluid dominated solutions 

Next, we consider the point (a6):(x, y, 7) = (0,0,0) 
which corresponds to a new type of fluid dominated so- 
lution. We don't consider is as one of the standard fixed 
points because of the apparently singular behaviour of 
the point in various limits. For example in Eqs. ([22]) - 
([24]) . x = is singular for q > — 1/2, y = is singular for 
q < and 7 = for q > 1/2. It complicates the analyis 
for this point but the procedure itself is similar to the 
case with (a5). 

Writing x = + Sx, y = + Sy, 7 = + ^7 and keeping 
only the leading terms, we obtain 



$y,N = 



Sx 
~2 

Sy 



3w r . 



%/3A 



(Sy) 



2q+2 



Sj,N = 



3(1 + w m ) - v^A 
3 + V3A 



(6xfi+ 1 (6i)i-i. 



(6y) 2 » (r (Sy) 2 



(95) 
(96) 



(Sx) 2 i- i (s^y-^ \p (Sx) 



(97) 



Following the previous case (a5), we introduce 9 = 
r and (3 = g^. Eqs. ([95])- (07]) then yield: 



9 m 3 



2 (l + Wm ) + ^-~X9 



-2q[ 1 + -) +- 
pj p 



& = 3-V3A0(l + /3), 

which has three fixed points given by 
(0,/3) = (O,O) , and 

a/3(1 + Wm) 



A 



2q[l + 
a/3(1 + io m ) 



.0 



and 



(98) 
(99) 

(100) 
(101) 



(2 9 -l) 1 



2 ^( 1 + i)-? 



1 + w„ 



(102) 



For the eigenvalues corresponding to Eq. (|100l) we ob- 
tain 3(1 + w m )/2 and 3, which means this fixed-point is 
unstable. 

Considering the fixed-point Eq. (|101[) first, the condi- 
tions for 9 to be a finite constant are, q > 1/2 or q < 0. 
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If < g < 1/2, then in the definition of 9, all the terms 
5x, 5y and 5j would appear in the numerator, with no 
apparent singularity present, implying the solution would 
not be of the form being discussed in this section. In ad- 
dition to this, as 5x, 5y, 5^ are all positive, 9 must be 
positive, which gives the following condition 



(103) 



2q 1 + - - - > : 
PJ P 



which is r < —2 with p > 0, or r > —2 with p < 0. 

Small perturbations around the fixed point Eq. (llOip . 
lead to the following eigenvalues for 58 and 5(3, 



"-(I + W m ) 



wp = 3- 



3(1 + W m ) 

2q{l 



p) p 



(104) 
(105) 



where from Eq. (|103|) . wp < is satisfied as long as 

(l + u» m )>2?(l + -)-->0, (106) 
P P 

which is —2 — p(l + w m ) < r with p > 0, or —2 — p(l + 
w-m) > r with p < 0. 

Having shown these new fixed points exist and have 
well defined behaviour through their combination in 
Eq. (fTUTj) . we perturb about them in Eqs. {22} -([2l|) to 
obtain the following eignenvalues for the perturbations 
in x, y and 7 



W-2 



2 Wm ) 

3, . f 1 

= — ( 1 + W m ) 1 ; — 

2 K m, [ 2?(l + £)- 

3. ,/p + r + 2\ 



(107) 



(108) 



w 3 = 3 1 



' 2 — ru7 r) 
r + 2 



(109) 



where we have used the relation q = —l/(p+r). Although 
W\, W2 and w 3 diverge for r = —2, from Eq. (|103p . this 
case is excluded. 

Notice that although we have shown eigenvalues for 
3 directions, only two of them are independent of the 9 
direction. This can be seen by the fact that wi, u>2 and 
W3 are related through: 

(1 - 2q) Wl + 2qw 2 + (i - q)w 3 = . (110) 

We can summarise the constraints on the parame- 
ters if the solutions approaching the point (a6) given by 
Eq. (jlOip are stable. From q > 1/2 or q < (or equiva- 
lently > p + r > —2 or p + r > 0) and W2 < 0, we have 



r < —2, which means — 2— p(l+w m ) < r < —2 withp > 
from Eqs. (|103|) and (|106j) . From tui < and w 3 < we 
have 2/r < w m < 0. As p > 0, > 1 > r/p > — 2/p for 
q > 1/2 and 1 + r/p > for q < 0. 

We now turn to consider the fixed-point given by (| 102|) . 
In this case only q > 1/2 is allowed as Sx and Sy have the 
same time dependence. Furthermore, as j3 is a nonzero 
finite constant, in addition to (|103p . another condition 
given by 



(111) 



{2q - 1) ( 1 + - ] - w m > , 



which is r < —2 — p{l + w m ) with p > or r > —2 — 
p(l + w rn ) with p < must be satisfied. In fact this 
extra condition is more restrictive than the one given by 
Eq. (TTU31. 

By considering small perturbations around this point, 
we obtain the following eigenvalues for 59 and 5/3, 



w 



Wp 



■-(1 + Wm) 

-3[(2g-l)(l 



w n 



2<z(i + 5) 



(112) 
(113) 



which arc both negative once the conditions given by 
Eqs. (IT03"|) and (jPTTI) are satified. 

By substituting f|102[) into Eq. we obtain the fol- 
lowing eignevalue for the perturbations in x 



w 1 



3(l + ^)(-l + 2q)(l+w m ) 
2[2g(l + 
3(p + r)(2g- l)(l + u; m ) 
2(r + 2) 



(114) 



where we have used the relation q — — l/(p + r). 

Notice that we have already obtained the stability for 
the 9 and (3 directions, the stability based on the eigen- 
values for the directions of x, y and 7 have the same 
information. Actually, the eigenvalues corresponding to 
the perturbation in y and 7 are given by 



w 2 = 



w 3 



3(1 + p(-l + 2 g )(l + «; m ) 

2[2q(l + £)-£] 
6(l+£)(l+™ m ) 



2[2g(l + 
3(p + r)(l + w m ) 
(r + 2) 

which are consistent with 

Wi=W 2 , 



(115) 
(116) 
(117) 

(118) 



(-2q + l)w 1 +2qW2 + (--q)w 3 = 0, (119) 

derived from the fact that 9 and (3 are constant. Although 
u>i, u>2 and w 3 diverge for r = — 2, from Eq. (I103p . this 
case is excluded. 
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We can summarise the constraints on the parame- 
ters if the solutions approaching the point (a6) given 
by Eq. (|102j) are stable. From q > 1/2 (or equivalently 
> p + r > —2) and wi,2,3 < 0, we have r < —2, which 
means —2 — p(l + w m ) > r with p > from Eq. If 
we interpret this as the constraint on w m , together with 
q > 1/2, w m < should be also satisfied. As p > 0, 
> 1 + r/p > -2/p for g > 1/2. 

It is worth mentioning that for the fluid dominated 



solution (a6) realised both through Eqs. (jlOll) and (|102l) . 
u> m < is required. Such a restriction implies that we 
can not be considering the usual matter or radiation fluid 
as the background for (a6) to be stable under, although 
it is consistent for example with that of a cosmological 
constant. 

We finish this subsection by providing a summary of 
the solutions (al) to (a6) and their stability criteria in 
Table U below. 





X 


V 


7 


il$ 


Valid q 


Stability 


(al) 





1 


1 


1 


q < -1/2 and q > 


marginally stable 


(a2) 





1 





1 


-1/2 < q < 1/2 


marginally stable 


(a3) 


1 








1 


q < 


w m > 0, r < -2, p > 0, £ > -1 












q> 1/2 


w m > 0, r < -2, p > 0, -1 - § < T - < -1 


(a4) 


x (0 < x < 1) 








x'i 


q < 


w m = 0, r < -2, p > 0, | > -1 












q> 1/2 


w m = 0, r < -2, p > 0, -1 - | < J < -1 


(a5) 








1 





<2 < 


g§ >w m >-l, p>2, e>-1 


(q6) 














q < 


> ™ m > £ , p > 0, -2 > r > -2 -p(l + w m ), ~ > -1 












«> 1/2 


0>w m (>§), p>0, -2>r, _1_S<E<_1 



TABLE II: This table summarises the validity and stability of the fixed-points al to ad. Notice that we have restricted ourselves 
to the regime A > and x > 0. For the condition for stability of (a6) for q > 1/2, w m > 2/r holds only for r > — 2— p(l + w m )- 



E. Late-time behaviour 

Now that we have analyzed the stability of the fixed- 
point solutions, we are in a position to discuss the late- 
time attractor structure of the autonomous system de- 
scribed by Eqs. (|2"2"1) - (|2~i)) . It is clear from the solutions 
presented in TABLE fl] and [TT] that the precise structure 
of the solutions to the dynamical system depends on the 
value of q. We have seen that for the particular cases of 
q = 1/2, q = 0, q = — 1/2 there are non-trivial "scaling 
solutions" where x, y and 7 are finite constants depend- 
ing on the model parameters A, p,, and the equation of 
state of the fluid w m - Therefore, in this subsection, first 
we will discuss the late-time attractor structure of the 
system with special emphasis on these three cases, and 
then we will move on to a discussion of the other cases. 



1. The case with q = 1/2 

TABLE HI shows that in the case of X+p, < (or equiv- 
alently p > 0), if w m > 0, the standard fixed-point (a4) 
(ultra-relativistic kinetic-potential scaling solutions) is 

stable for w m > — 1 + X(^/X 2 — 12 — A)/6. It is also stable 
for w m < if A < \/3(l + ui m )/y/—w m , which also corre- 
sponds to the requirement that p > with p+r = —2. On 
the other hand, if A > y3(l + w m )/ \J—W m with w m < 



(or equivalently p > v / 3cw(l + w m ) / V '—w m > 0), then 
the standard fixed-point (a5) (ultra-relativistic kinetic- 
potential-fluid scaling solution) is the stable solution. 
Clearly, these two fixed-point solutions are candidates for 
the late-time attractor solution to the system as long as 
the model parameters satisfy the conditions mentioned 
above. On the other hand, TABLE [Til shows that for 
q = 1/2, the fixed-point associated with (al) (standard 
potential dominated solutions) is marginally stable. As 
(al) has a zero eigenvalue, to specify which of the three 
possibilities (<z4), (a5) or (al) is likely to win out will 
require that we go to higher order, although it is likely 
to be (o4) or (a5) as one of these is definitely stable when 
the conditions A + p < and 9=1/2 are satisfied. It 
is worth mentioning that for the cosmologically sensible 
cases, we expect w m > 0. Then, only (a4) can be the 
late-time attractor. 

Alternatively, for A + fi > (or equivalently p < 0), 
as both points (a4) and (a5) have positive eigenvalues, 
we can expect (al) to be the late-time attractor. These 
findings have been confirmed numerically as can be seen 
in Fig. [Q 



As we have already mentioned, (a4) and (a5) are 
new solutions found only in the case of the DBI field. 
(a4) corresponds to the ultrarelativistic version of the 
usual canonical power-law inflationary solution. It is well 
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FIG. 1: Demonstration of the attractive property of (ad) in 
the regime A + ji > and q — 1/2. Note that for all the 
initial values of x, y and 7, the solutions approach that of al 
where (x,y,7) = (0, 1, 1). Note also that the physical region 
is conned to x 2 + y 2 < 1 in the three-dimensional phase space. 
Although here we have chosen q — 1/2, r/p — —0.5, A = 10 
(or p — —4, r = 2, (av) 1 / 2 = 2/5), w m = 0, the same results 
are obtained for q = 1/2 with r/p > — 1. 




I 

FIG. 2: Late-time attractor structure in the models where 
the potential and the brane tension are power-law functions 
of the DBI field with q — 1/2 and w m > 0. We show the 
late-time attractor solution for a given model parameter in A 
- fx space with A > 0. In the case with w m < 0, the region 
where (a5) is the late-time attractor appears for A + fx < 
and A > V3(l + w m ) / 'y/—w m - It is worth noting that A and 
fx are related with the parameters of the potential and brane 
tention as A = \p\/ '-Jav and fx = — (2 + p)\p\/ (p^Jov). 



known that power-law inflation is strongly constrained by 
the observation of the spectral index of primordial pertur- 
bations (for example, see [HJ ). Although it is difficult to 
have DBI driven power-law inflation, we can in principle 
use it to explain the current acceleration of the universe. 
On the other hand, although (a5) is also unique to the 
DBI field, because the solution requires w m < to be 
satisfied, it is not so realistic given that we are usually 
considering w m to be either matter or radiation. 

However, since (al) corresponds to the usual acceler- 
ating solution where the scalar field rolls slowly down its 
potential and is not unique to the DBI field, this solution 
can in principle be used to explain both an early stage of 
inflation and the present dark energy dominated period 
of acceleration. 

In Fig. [3J we summarise the late-time attractor struc- 
ture with q — 1/2 and w m > 0. 



2. The case with q = 

Here, we consider the late-time attractor structure in 
the models where the potential and the brane tension are 
power-law functions of the DBI field with q = 0. 

As is shown in TABLE U in the case with A + p, > 
(or equivalently 1 + r/p > 0), if A < y3(l + w m ) (or 
equivalently \p\ < y/3(l + w m )), the fixed-point (62) (the 
standard kinetic-potential scaling solution) is stable. If 
A > y/3(l + w m ) (or equivalently \p\ > + w m )), 

the fixed-point (63) (standard kinetic-potential-fluid scal- 
ing solutions) is stable. Furthermore, if the condi- 
tions p, < — v6 and w m > 3/(/i 2 — 3) (or equivalently 
w m > 3/(r 2 — 3)) are satisfied, (cl) (the relativistic ki- 
netic energy dominated solution) is stable. Therefore, 
these fixed-point solutions are clearly candidates for the 
late-time attractor behaviour, provided the model pa- 
rameters satisfy the conditions mentioned above. We 
also find for the region in the parameter space satisfy- 
ing both stability conditions for (63) and (cl), which is 
the late-time attractor depends on the initial values of 
x, y and 7. On the other hand, TABLE UH shows that 
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I 

FIG. 3: Demonstration of the attractive property of (a2) in 
the regime \ + p, < and q = 0. Note that for all the ini- 
tial values of x, y and j, the solutions approach that of (q2) 
where (a;, y, 7) = (0,1,0). Note also that the physical re- 
gion is conned to x 2 + y 2 < 1 in the three-dimensional phase 
space. Although here we have chosen q = 0, r/p — —1.5, 
A = 10, Wm = 0, the same results are obtained for q — with 
r/p < — 1. Notice also that we chose the directions of the axis 
differently from Fig. [1] 

I 



for <7 = 0, the fixed-point (a2) (ultrarelativistic potential 
dominated solutions) is marginally stable. 

Using arguments similar to those applied for the case 
with q = 1/2, we expect that for A + p, > (or equiv- 
alently 1 + r/p > 0), one of (62), (63), (cl) will be the 
late-time attractor solution of this system as long as their 
stability conditions are satisfied. This is because we have 
seen that the stability properties of these three points are 
given by three negative eigenvalues which is a stronger 
condition than for (a2). 

On the other hand for A + fi < (or equivalently 
1 + r/p < 0), the points (62), (63) and (cl) have positive 
eigenvalues, so we expect (a2) will become the late-time 
attractor. This has once again been confirmed numeri- 
cally in Fig. |3l 



Since (62) and (63) correspond to the well known 
power-law inflationary and scaling solutions respectively, 
the cosmology based on these solutions has already been 
well studied [53T - [57l | ■ As in the ultra-relativistic case, 
the power-law inflationary solution (62) can explain the 
present acceleration, while the scaling solution (63) can 
play a very important role in classifying the late-time 
attractor structure of the system. 

As we have mentioned earlier, (cl) is a new solution 
specific to the DBI field, but we do not consider it further 
because it can neither explain the current acceleration of 
the Universe, nor can it accommodate a fluid. 

(a.2) is also a solution specific to the DBI field and is in 
fact a very interesting new type of inflationary solution. 
These models are however tightly constrained. They have 
been shown to give too large a degree of primordial non- 
Gaussianity 0-11]. 

In Fig. 21 we summarise the late-time attractor struc- 
ture with q = and w m > 0. 



3. The case with q = - 1/2 

Here, unlike the previous cases with q — 1/2 and q = 0, 
when q = —1/2, the stable solutions with three negative 
eigenvalues for X+p > (or equivalently p > 0), are given 
by the four fixed-points ((a3) - (a6)) shown in TABLEQl] 

(a3) (ultrarelativistic kinetic dominated solutions) is 
stable for w m > with p > 4, while (a4) (Ultrarelativis- 
tic kinetic-fluid scaling solutions) is stable for w m = 
with p > 4. (a5) (Standard fluid dominated solutions) 
is stable for (p — 6)/(p + 2) > w m > —1, with p > 2, 
and (a6) (Ultrarelativistic fluid dominated solutions) is 
stable for > w m > 2/(2 - p), with A/(-w m ) > p > 4. 

These four fixed-points are candidates to be the 



late-time attractor solutions to the system as long as 
the model parameters satisfy the conditions mentioned 
above. When a region of parameter space allows more 
than one stable fixed-point, which one of them becomes 
the late-time attractor depends on the initial values of x, 
y and 7. On the other hand, from TABLE HI we see that 
in this case the fixed-point (c3) (the relativistic potential 
dominated solution) is marginally stable. Therefore (c3) 
is also a candidate late-time attractor solution. However, 
as we saw before, when a marginally stable solution also 
exists, we expect that for A + p, > (or equivalently 
p > 0), one of the four fixed points (a3)-(a6) will be 
the late-time attractor solution of the system in the re- 
gion of the phase space when their stability conditions 
are satisfied. 
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V3(1+Wra) 




FIG. 4: Late-time attractor structure in the models where the potential and brane tension are power-law functions of the DBI 
field with q — and w m > 0. We show which fixed-point solution is the late-time attractor solution for a given set of model 
parameters in X-jl space with A > 0. In the case with w m < 0, the region where (cl) is the late-time attractor disappears. This 
figure can be also applied to the models where the potential and brane tension are exponential function of the DBI field with 
q — through the replacements A(= |p|) —¥ A, /x(= r\p\/p) — > /x. For this class of models, it is possible to consider the case with 
fi = —A where (c4) or (c5) can be the late-time attractor. Although it is complicated to write down, there is a critical value of 
A depending on w m , a and v, below which (c4) is the late-time attractor and above which (c5) is the late-time attractor. 



Turning our attention instead to consider the regime 
A + Jh < (or equivalently p < 0), we see that the fixed- 
points (a3) - (a6) now have positive eigenvalues, and so 
we expect (c3) to become the late-time attractor, a result 
that we have confirmed numerically in Fig. [5J 



Cosmology based on the existence of the solutions (a3) 
and (a4) are peculiar to the DBI field in the sense that 
the field behaves like dust even though the kinetic term 
completely dominates the potential term. However, al- 
though interesting in its own right, the solutions do not 
appear to be useful cosmologically. Solutions (a5) or 
(a6) both imply that the Universe is completely domi- 
nated by the background fluid with the DBI field playing 
a negligible role in its evolution. 

On the other hand, the solution (c3) appears to have 
some interesting features specific to the DBI field. In 
this case, 7 or equivalently the sound speed c s can be 
a constant ranging between and 1. In such a case, 
if this solution is realised in the very early Universe, it 
opens up the possibility of a large but still allowed de- 
gree of non-Gaussianity being obtained. Furthermore, 
since the current prediction of non-Gaussianity from this 
type of inflation model is based on the assumption that 
c s is constant, this new solution serves as a good con- 
crete background about which to consider cosmological 
perturbations. 



4- The case with general q 7^ 0, 1/2 or —1/2 

To complete the classification of the late-time attractor 
structure in the models where the potential and the brane 
tension are power-law functions of the DBI field, we now 
consider the case where q is different from 1/2, or —1/2. 
For three cases we have earlier seen that the question of 
whether marginally stable solutions can be the late-time 
attractor or not depends upon the stability of the other 
allowed fixed-points. In particular those fixed points with 
three negative eigenvalues, means that these marginally 
stable fixed-point solutions can not be the late-time at- 
tractor as their stability is weaker than that of the fixed- 
points. However, if all other fixed-points have positive 
eigenvalues, then the marginally stable fixed-point solu- 
tions can turn out to be the late-time attractor. 

Significantly, since this criteria also holds for the cases 
considered here, we will apply it without showing the 
numerical plots except the special case with < q < 1/2. 

For q > 1/2, for A+/2 < (or equivalently p > 0), there 
are four stable fixed-point solutions (a3) — (a6) . In addi- 
tion to these, the fixed-point solution (evl) is marginally 
stable. 

As in the previous cases, for X + jl < 0, one of the fixed- 
point solutions (a3) — (a6) can be the late-time attractor 
depending on the values of parameters p, r, w m or the 
initial values of x, y, 7, while for X+fl > (or equivalently 
p < 0), (cvl) is the late-time attractor. 

For < q < 1/2, there are two marginally stable fixed- 
point solutions (al) and (a2). As there are no other 
fixed-points whose stability around the points are charac- 
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FIG. 5: Demonstration of the attractive property of (c3) 
in the regime A + ji < and q = —1/2. For every ini- 
tial values of x, y and 7, the solution approaches (a;, y, 7) = 

(0,1, \/3/-\/A 2 + 3). Notice that the physical region is con- 
fined to x 2 + y 2 < 1 in the three-dimensional phase space. Al- 
though we choose q — —1/2, r/p — —1.5, A = 10 (or p — — 4, 
r = 2, ((j!/) 1//2 = 2/5) and w m = in this example, the same 
results are obtained for q = —1/2 with — 1 > r/p. Notice that 

for A = 10, VH/y/x 2 +3 ~ 0.17 as expected. 

I 



terised by three negative eigenvalues, determining which 
of the two will be the late-time attractor is non-trivial. 
However, numerically we find for A + jl > (or equiva- 
lently p < 0), (al) is the late-time attractor whereas for 
A +/z < (or equivalently p > 0), it is (a2). (See Figs. [5] 
andO 

In the case where —1/2 < g < 0, when A + jl > 
(or equivalently p > 0), there are four stable fixed-point 
solutions (a3) — (a6), with the fixed-point solution (a2) 
being marginally stable. One of (a3) — (a6) can be the 
late-time attractor depending on the values of the param- 
eters p, r, w rn or the initial values of x, y, 7, while for 
A + /i < (or equivalently p < 0), (a2) is the late-time 
attractor. 



When q < —1/2, for A + (I > (or equivalently p > 0), 
the story is same as in the case with > q > —1/2, there 
are four stable fixed-point solutions (a3) — (a6). But 
in this case, an additional marginally stable solution is 
(al). Once again, one of (a3) — (a6) can be the late-time 
attractor solutions depending on the values of parameters 
p, r, w m or the initial values of x, y, 7, while for X+p. < 
1 1 




1 



FIG. 6: Demonstration of the attractive property of (al) in 
the regime A + jl > 0. Note that for all the initial values of x, 
y and 7, the solutions approach that of al where (x, y, 7) = 
(0, 1, 1). Note also that the physical region is confied to x 2 + 
y 2 < 1 in the three-dimensional phase space. Although here 
we have chosen q — 1/4, r/p — —0.5, A = 10 (or p = —8, 
r = 4, (au) 1 ^ 4 = 4/5), w m = 0, the same results are obtained 
for q > with r/p > —1 and —1/2 > q with r/p < —1. 
I 



(or equivalently p < 0), (al) is the late-time attractor. 

We summarise the possible late-time attractor solu- 
tions for a given set of q, p, r based on the discussion in 
this subsection in TABLE [TTT1 



IV. MODELS WITH EXPONENTIAL the DBI field, 

POTENTIAL AND BRANE TENSION 

In this section we consider the models where both the 
potential and brane tension are exponential functions of 
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1+ -> 


1 + - < 




/ 1 \ 

(al) 


(aS) or (a4) or (ao) or (ab) 




(al) 


(a4) or (a5) 


y > q > 


(al) 


(q2) 


g = 


(62) or (63) or (cl) 


(a2) 


0>g>-± 


(ct3) or (a4) or (a5) or (q6) 


(a2) 




(a3) or (q4) or (a?5) or (q6) 


(c3) 


-5 > 1 


(a3) or (q4) or (a?5) or (a6) 


(al) 



TABLE III: This table summarises which fixed-points can be the late-time attractors for the models where the potential and 
brane-tension are power-law functions of the DBI field for a given q, p and r. "or" means there may be more than one possible 
late-time attractor and which is realised depends on w m and the initial value of x, y, 7, which has been described in TABLE U 
and TABLE HQ in detail. 




A. Autonomous System 

In this case, from the discussions in Sec. HH in terms 
of constants A and n given by Eq. ([120]) . Eqs. (TTSJ) - (fT5]) 
constitute an autonomous system. 

These equations correspond exactly to Eqs. (|22|) - (f24)) 
with g = under the identification A — > A(= \p\) and 
/j — > fi{= r\p\/p). It is as expected because the exponen- 
tial function can be regarded as the power law function 
cj) p with p — > 00. This means that all the fixed points 
associated with q — are found in the case of the expo- 
nential functions as well. However, as we shall now show, 
when A + [i = 0, there is an additional fixed point found 
in the exponential case. As in the previous section, with- 
out loss of generality, we restrict our discussion to the 
case A > and x > 0. 



B. Fixed-point Solutions 



FIG. 7: Figure showing the attractive nature of the fixed 
point (a2) in the regime A + fi < and q ~ \. For all initial 
values of x, y and 7, solutions approach the point (q2) where 
(x,y,y) = (0, 1,0). Notice that physical region is confined to 
x 2 + y 2 < 1 in the three-dimensional phase space. Although 
we have chosen q — 1/4, r/p = —1.5, A = 10 (or p = 8, 
r — —12, (a-;/) 1 / 4 = 4/5), w m = in this example, the same 
results are obtained for 1/2 > q > —1/2 with —1 > r/p. 
Notice also that we chose the directions of the axis differently 
from Fig. Q] 



r 



V{4>) = ae 



-X(j> 



(120) 



where a and v are constants. We will see that there is a 
direct link with the q = limit for the scaling solutions 
(b2), Eq. ((3HJ> and(b3), Eq. (J35J) obtained in the previous 
section. 



From our discussion in the power-law case for the po- 
tential and brane tention with q = 0, the fixed-point solu- 
tions obtained there are also fixed point solutions for the 
exponential case. This means that the eight fixed-point 
solutions (al), (a2), (61) -(63), (cl), (c2),(a2) obtained 
in Sec. Mil are also fixed-point solutions with the replace- 
ment A(= \p\) — > A and Jx{= r\p\/p) — > /i. 

The additional solution not present in the power law 
case follows when 



jj, = —A . 



(121) 



As we will now show a class of new solutions then exist 
for < 7 < 1 and non-zero x and y . 

Strictly speaking, since in this case V(<f>) and f(<j>) are 
related through the constraint fV = ov, this is a two- 
dimensional system as opposed to a three-dimensional 
one. Actually 7 can be completely specified once x and 
y are given, through the relation 



7 



y 



3aux 2 + y 2 



(122) 



However, for simplicity we continue to solve explicitly the 
evolution equation of 7 to obtain the fixed-points. Of 
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course, from Eq. (|122[) . 7 becomes constant when x and 
y become constant, which provides a consistency check 
for the method. In fact the constraint Eq. (|122j) . leads 
to some more interesting constraints on the DBI field. 
Substituting it into Eq. (fTTj) and defining r = y 2 /x 2 , the 
equation of state of the DBI field is given by 



wj,(r) 



r(l — 3av — r) 
(3av + r)(l + r) 



with < r , (123) 



which clearly asymptotes to as r — > 0, while —1 as 
r — > 00. In fact the function decreases monotonically 
and is constrained as 



1 < W0 < . 



(124) 



To find the fixed-points in the system, substituting 
Eq. (|12ip into Eq. (|15[) , and requiring 7 remain constant 
implies, 



1 



1 



a/3(1 + 7) , 
V7A 



(125) 



Substituting Eqs. ([T2T1) and (|T2"5l) into Eq. (O, leads 
to the following two constant solutions for x: 



V7A 



\/3(H-7) 



or x 



\/3(l + W m ) 

^7(1 + 7) A' 



(126) 



The first of these solutions leads to x 2 + y 2 = 1, which 
from Eqs. (|9|) and (|10[) implies that during scaling Q m = 
and f2^ = 1. It corresponds to the relativistic kinetic- 
potential scaling solutions given by 



(c4) (x,y,7) 



V7A 



7A 2 



v / 3(r+^) 



3(1 + 7) 



■ 7 



with0< 7 <l, A < J ^t^ , [i 



-A. 
(127) 



From Eqn. (1111) it implies w§ — 1 which means 

that the solution is accelerating if A 2 7 < 2. Since for this 
solution, the DBI field completely dominates the fluid, 
this solution exists irrespective of the value of w m . 

Although somewhat complicated, the constraint 
Eq. (|122l) . allows us to determine the actual values for 
7, x and y for this fixed-point solution in terms of the 
parameters A, er, v. For example, for av > 2/3 they are 
given by 



7 



A 2 + ^36 + 12(-1 + iav)\ 2 + A 4 ' 

2A 2 

6 + A 2 + v /36 + 12(-l + 3cr^)A 2 + A 4 



1 1 2 

y = l-z , 



(128) 

,(129) 
(130) 



while for av < 2/3, a related but slightly different ex- 
pression holds. 

The second solution in Eq. (|126j) corresponds to the 
relativistic kinetic-potential-fluid scaling solution and is 
characterized by 

(c5) (x,y,7) 



V3(l + W m ) / 3(1 + W m )(7 - W m ) ~ 

VW+T)*1 ^7(1 + 7) ' 7 

with < 7 < 1, w m < 0, 

X>JW + Wm) , M = -A. 



(131) 



From Eqns. (fT0|) and ([TT|) we see that it leads to = 
3( - 1 ^"'- > and = w m . In other words the DBI field 
equation of state tracks that of the background matter. 
Since the property of this solution is completely affected 
by the nature of the fluid, the existence condition of this 
solution depends on the value of w m - It is worth noting 
that when we take into account the constraint Eq. (|122l) , 
then from Eq. (|124[) . < 0. Therefore, although it 
cannot be seen from Eq. (I131[) explicitly, this solution 
exists only for the case with w m < 0. 

As in the case of (c4), we can solve explicitly for 7, 
x and y in terms of A, a, v, w m by making use of the 
constraint given by Eq. (|122[) . For example, for the case 
with av > 1/3, 



-2w r< 



(3av - 1)(1 + w m ) 
1 



4m m 1 1 



y = 



' 1 ' (l-3<7I/) a (l + U) m ) 

18av(l + w m ) 
{3av-l){l + w m )X 2 



(132) 



1 



x 2 = 



(133) 
(134) 



1 + w„ 



7 - w„ 



-y , 



while for av < 2/3, a slightly different expression holds. 

In summary we see that there are ten fixed-point so- 
lutions for the case where both the potential and brane 
tension are exponential functions of the DBI field. In 
TABLE IIVI we summarise the two solutions (c4) and 
(c5) which appear when the particular condition A = — [i 
is satisfied. The other solutions correspond to the case 
q = and are given in TABLE HI 
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a/3(1 + 7) 



Existence 



Stability 



(c4) 



7A 2 
3(1+7) 



< 7 < 1,A < 



3(1+7) 



/i = — A, VW„ 



A < 



3(1+1U,„) 

7 



0<7< l,w m <0,A> ^HiH£ 2 a ) ^ = _A 



a(> 0) 



(c5) 



V3(l+-uj m ) 

y 7 (i+7)A 



3(l + m m )(7-u) m ) 
A 2 7(l + 7) 



3(l + u; m ) 
7T 2 



TABLE IV: Summary of the extra fixed-points and their stability in the models where the potential and brane tension are 
exponential functions of the DBI field satisfying A = — fj,. It can be shown that for (c5), if the existence conditions are satisfied, 
then the stability conditions are automatically satisfied. Although we keep 7, from the constraint given by Eq. (|122|l . it is 
specified in terms of A, a, v, w m through Eq. (|128|1 for (c4) and Eq. (|132p for (c5). In addition to these two, the solutions 
(al), (o2), (61) — (63), (cl), (c2), (a2) obtained in Sec. IIII Bl are also fixed-point solutions with the replacement A(= \p\) — > A 
and p,(— r\p\/p) — > [i. Notice that we have restricted ourselves to A > and x > 0. 



C. Stability Analysis 



potential-fluid scaling solutions) we obtain: 



We turn now to study the stability of the fixed-point 
solutions obtained in the previous subsection. Given 
that eight of the fixed-point solutions ((al), (a2), (61) — 
(63), (cl), (c2), (a2)) were also obtained and their stabil- 
ity determined in Sec. IIIII we do not repeat that anal- 
ysis here. It is basicially the same as in that case 
except we make the replacement A(= \p\) —¥ A and 
//(= r\p\/p) — > //. Here, we concentrate on the two ex- 
tra solutions (c4) and (c5) in Eqs. (|127[) and (|13ip . For 
the fixed-point (c4) (relativistic kinetic-potential scaling 
solutions), having perturbed about the solution to linear 
order, we obtain for the corresponding eigenvalues: 



101,2 = |(1 - w m ) 



A^-6 



w 2 = A 7 — 3(1 + w m ) 

W3 = . 



(135) 



W3 is zero because this corresponds to the 7 direction and 
at this fixed point, the RHS of Eq. (| 15[) becomes 0. As 
mentioned earlier, this is to be expected because there 
is a constraint on 7 and only two of the three degrees of 
freedom are actually independent for A = — ji. Therefore 
the existence of two negative eigenvlues will guarantee 
the stability of this fixed-point solution. 

/B7, W\ is always negative. 



Since A < x /3(l+7)/7 < 
Therefore, solution (c5) is stable if W2 < 0, that is, A < 
a/3(1 + Wm)/^- Recall that the condition for inflation in 
this case is A 2 7 < 2. 

It is worth mentioning that by adopting the constraint 
given by Eq. (|122[) . actually 7 is completely specified in 
terms of A, cr, v. For example if av > 2/3, it can be 
shown that this solution is stable if 

[6(3 + (3av - l)A 2 )]w m + 18 - 6(1 - 3av)X 2 + A 4 
> v /36~12(l-3cr^)A 2 + A 4 A 2 . 

(136) 

Similarly, for the fixed-point (c5) (relativistic kinetic- 




8(1 + w m )(w m - 7 2 )(3 + 3iu m - A 2 7) 



A 2 (-l + w m ) 2 7 



(137) 



The same reasoning as applied to (c4) suggests that 
even though w 3 = 0, the solution (c5) is stable if the 
real parts of both w\ and W2 are negative. This requires 
w m < 7 2 , which is automatically satisfied if the existence 
condition w m < is satisfied, although it would imply an 
unusual form of background matter being considered. 

As in the case of (c4), by adopting the constraint given 
by Eq. (|122l) . actually 7 is completely specified in terms 
of A, cr, v, w m - It can be shown that if av > 1/3, the 
existence condition A > y3(l + w m )/Vl which is the 
most stringent condition for this solution to be stable 
can be expressed as 

[6(3 + (3cn/ - l)A 2 )]w m + 18 - 6(1 - Zav)\ 2 + A 4 



< ^36 - 12(1 - 3cr i /)A 2 + A 4 A 2 . 



(138) 



Clearly, Eqs. (fl"3"o) and (ITSS]) show that for av > 2/3, 
there is no overlap in the parameter space between where 
(c4) and (c5) are the late-time attractors. In TABLE HVl 
we summarise the stability conditions of these two fixed- 
point solutions. 

There has been related work looking at the fixed-point 
solutions (c4) and (c5). Actually, (c4) was previously 
obtained in }50l . [52| . while the existence of the scaling 
solution (c5) was pointed out by Martin and Yamaguchi 
in [4911 . The authors' stability analysis followed that of 
,54], which although providing a proof that the fixed- 
point solution is attractive, it did not explain how the the 
solution could be realised from all initial values. In our 
approach we have gone into more details, following the 
analysis of 57]. In particular by making the phase space 
compact, we have been able to establish all the fixed- 
points in the compact phase space. As we will now see it 
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then becomes possible to discuss the late-time attractor 
structure, something we now turn our attention to in the 
following subsection. The solution (c5) has recently led 
to a number of papers investigating its cosmology [49-52] 

D. Late-time behaviour 

In Sec. IIII1 Fig. [U shows the late-time attractor struc- 
ture with q — and w m > 0. From the discussion at the 
beginning of this section, except the special case with \x = 
— A, Fig. S] shows also the late-time attractor structure in 
the models where the potential and brane tension are ex- 
ponential functions of the DBI field with the replacement 
-H= \p\) ^ an d A( = r \p\/p) M- I 11 the case with 
fi > —A, for A < y3(l + w m ), the fixed-point (62) (stan- 
dard kinetic-potential scaling solutions) is the late-time 
attractor, whereas for A > + w m ), the fixed point 

(63) (standard kinetic-potential-fluid scaling solutions) is 
a late-time attractor. For /i < — \/3(l + w m ) / yjw m with 
w m > 0, in addition to (63), there is a possibility that 
the fixed-point (cl) (relativistic kinetic dominated solu- 
tions) is also a late-time attractor. Because, they are all 
stable locally, which of these solutions actually wins out 
depends on the initial value of x, y and 7. 

On the other hand, from the discussions in Sec. IIII1 
when fi < —A, the fixed-point solution (a2) (ultra- 
rclativistic potential dominated one) is the late-time at- 
tractor. For the case \x = —A, as we have just seen, the 
additional fixed-points solutions (c4) and (c5) are also 
possible late-time attractors, although the conditions un- 
der which they become attractors do not overlap in the 
available parameter space, so they do not compete with 
one another for overall stability. 

Of course when considering the overall stability for the 
case of /i = — A, we need to also include a discussion of the 
fixed-point solutions which can be late-time attractors for 
the general case /1 7^ A. Since (cl) and (a2) do not satisfy 
the constraint given by Eq. (|122p . these two fixed-points 
can not be the late-time attractor for /1 = —A. It can be 
shown that (c4) and (c5) reduce to (62) and (63), in the 
limit of 7 — > 1. 

It follows that only (c4) and (c5) are the late-time at- 
tractors for fx = —A and once we specify the values of the 
parameters, A, <r, v and w m we can judge which of these 
two fixed-points will be the late-time attractor. 

E. Power-law models with p + r = 

Here, for completeness, it is appropriate to mention 
what happens when the models V((f>) and /(</>) are given 
by 

V (</>) = *\<t>\ p , /(0)=^|- p . (139) 

From (11391) . as in the exponential potential case with 
/j, = — A discussed in section llVBl there is a constraint 



given by Eq. (|122[) . implying that only two of the three 
variables arc independent. This in turn implies that we 
cannot make use of the degree of freedom correspond- 
ing to f((j>) to construct an autonomous system to solve 
Eqs. (H~3]) - (TT6|) as they stand. So if there is a constraint 
of the form f(<f>)V(<f>) = const., the only case we can ob- 
tain an autonomous system is when V ~ exp[— X(f>] and 
/ ~ exp[A0] (or equivalently p — > 00). 

We can of course still make some progress. For example 
for the case of a canonical scalar field, not described as 
an autonomous system, the late-time behaviour with a 
power- law potential has been determined in [83 85] . We 
adopt a similar procedure here for the case described by 
Eq. (|139[) . by regarding these cases as limits of the models 
where the potential and brane tension are exponential 
functions of the DBI field satisfying fj, = — A. The point 
is that in terms of A and /z, the evolution equations for 
x, y, 7 are still given by Eqs. ([T3]). (fT4|) and ([15]). The 
difference from the cases with the exponential potential 
and brane tension is that A and /1 are not constant for 
the power-law cases but are given by 

X = -fi = -ep^- (140) 

m 

As in the other cases, we can generally restrict A > 
which is equivalent to considering only (f> < for p > 
and (j) > for p < 0. Therefore, in the late-time limit, 
<p — > for p > while 4> — > 00 for p < 0. This means 
that the late-time asymptotic value of A is A — > +00 for 
p > 0, while A -> -t-0 for p < 0. 

First, let us consider the case with p < (A — > 0). 
In this case, from Table ITVl and Fig. SI the only possible 
late-time attractor solution is (c4) (relativistic kinetic- 
potential scaling solutions). In the limit A — > 0, from 
Eqs. (TT2"g]) - ([T30]) , it asymptotes to (x,y,j) = (0,1,1), 
that is, standard potential dominated solutions. There- 
fore, the DBI-field behaves like a cosmological constant 
at late-time. 

Next, let us consider the cases with p > (A — >• 00 ). 
In these cases, since A is very large, from Table IIVI and 
Fig. 21 the only possible late-time attractor solution is 
(c5) (relativistic kinetic-potential-fluid scaling solutions). 

In the limit A — »• 00, from Eqs. (|132p - (|134[) . it asymp- 
totes to [x,y, 7) = (0,0,7), that is fluid dominated so- 
lution and the DBI field becomes irrelevant compared to 
the fluid. It is worth noting that from Eq. (|132p 7 is 
completely specified once er, v, w m are fixed. 



V. SUMMARY 

Successful models of inflation arising within string the- 
ory from the DBI action have generated a great deal of 
interest recently. With their non-canonical kinetic terms, 
non-trivial potentials and brane tensions, they have led 
to a number of fascinating results including the predic- 
tion of distinctive non-Gaussian fluctuations in the CMB. 
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Although most models investigated to date have had spe- 
cific functional forms for the potential and brane tension, 
in this paper we have decided to broaden the class of 
models being discussed and so have analysed the dynam- 
ics associated with more general forms for these poten- 
tial and brane tension functions, including in the analysis 
the presence of a background perfect fluid in a flat FRW 
universe. Following the approach developed in [57| . we 
have introduced a suitable set of dynamical variables x, 
y and 7 in Eq. ((HJ which has allowed us to determine 
the phase-space portrait of the system. In particular, we 
have established the late time behaviour of these systems, 
demonstrating where appropriate the attractor nature of 
the solutions. 

In Sec. IIII1 we have considered the models where the 
potential and brane tension are given by power-law func- 
tions of the DBI field (V(4>) = a\<j)\ p , /(<£) = v\<j)\ r ). The 
standard fixed-point solutions of this system are sum- 
marised in TABLE HI where we see that the late-time at- 
tractor nature of the solutions depends ong= — l/(p+r). 
The interesting cases of scaling where the ratio of the ki- 
netic to potential energies of the DBI field is a constant 
are found to exist only for q = l/2 ((o4), (a5)) and q = 
((62), (63)). This is because if we require x and y to be 
nonzero constants, there are only two possibilities, that 
is, 7 = (leading to scaling with q = 1/2) and 7=1 
(leading to scaling with q = 0). These scaling solutions 
are then shown to be stable for certain regions of the 
parameter space. In addition to these, we have also ex- 
plicitly demonstrated the existence and stability of an 
interesting inflationary solution (c3) specific to the DBI 
field with constant 7 which differs from 7 = 1 and 7 = 0. 

The DBI system is rich. For example the evolution 
equations (|2"2" j) -(|M )l can appear singular when some of x, 
y or 7 either tend to zero or unity, which one is singular 
depends on the value of q. On the face of it, the equa- 
tions appear to be ill-defined, but in practice it turns 
out that these points can actually be late-time attrac- 
tor solutions for the system. In section MI Dl we obtain 
these fixed points and determine their stability. These 
are summarised in Table [TTJ 

Having established all the fixed-points and their stabil- 
ity for a given set of parameters, we have gone on to de- 
termine which of these solutions will be the late-time at- 
tractor. Particular care is required when considering the 
cases where the eigenvalues associated with the pertur- 
bations vanish, implying that the solution is marginally 
stable. Whether these are the late time attractors for the 
system depends upon the stability of the other allowed 
fixed-points. In particular those fixed points with three 
negative eigenvalues, means that these marginally stable 
fixed-point solutions can not be the late-time attractor 
as their stability is weaker than that of the fixed-points. 
However, if all other fixed-points have positive eigenval- 
ues, then the marginally stable fixed-point solution can 
turn out to be the late-time attractor. We have sum- 
marised the possible late-time attractor solutions for a 
given set of q, p, r in Table III. 



In Sec. IIV1 we have considered the models where the 
potential and brane tension are exponential functions of 
the DBI field (V((f>) = Ce"^, /(<£) = De~^). This 
system has a similar dynamical structure to that of the 
power-law models with q — 0. However, there are ad- 
ditional scaling solutions present in the exponential case 
because we can construct an autonomous system for the 
case with A = — fj,, leading to solutions where 7 is a con- 
stant between and 1 ((c4),(c5)). In this case, we find 
that once we specify the values of the parameters, A, 
a, v and w m we can judge which of these two fixed- 
points will be the late-time attractor. The stability of 
these two fixed-points are summarised in Table IV. We 
have also shown that the special case (V(4>) = <i\(f>\ p , 
/(</>) = ^|</>|~ p ) can be discussed in terms of the limiting 
behaviour A — ¥ or A — > 00 in the models with an expo- 
nential potential and brane tension satisfying A = — fi. 

There is an overlap between elements of this work and 
other published material. For example, we are able to 
reproduce a number of results obtained earlier in [g9| 
where the authors considered the case with a massive 
potential and AdS throat. It turns out to be a special 
case of q = —1/2 discussed in Sec. Ill in this paper. 

The existence of the scaling solution (c5) was origi- 
nally pointed out by Martin and Yamaguchi in [49[ . The 
authors stability analysis followed that of J54[, which al- 
though providing a proof that the fixed-point solution is 
attractive, it did not explain how the the solution could 
be realised from all initial values. In our approach we 
have gone into more details, following the analysis of 
[57| . In particular by making the phase space compact, 
we have been able to establish all the fixed-points in the 
compact phase space, allowing us to then properly dis- 
cuss the late-time attractor structure of the solution. 

Of most interest to us though is the existence and sta- 
bility of a class of cosmologically relevant solutions. We 
have found that a fixed point solution (c3) (relativistic 
potential dominated solutions) where 7 is a constant sat- 
isfying < \/3/\/A 2 + 3 < 1 can be the late-time attrac- 
tor for q = — 1/2 with r/p < —1 (see also [EH). Although 
it can not be applied to the AdS throat (r = —4) as 
r/p > —1, we believe this solution is very interesting and 
important. In calculations of primordial perturbations of 
the DBI inflation models, properly incorporating the time 
de pen dence of the sound speed c s is a complicated issue 
(H |86|. and in fact it is usually assumed to be constant. 
(For a recent approach to relaxing this assumption, see 
87].) Since in our case c s — 7 is constant for the fixed- 
point solution (c3), this serves as a good background to 
use when we consider cosmological perturbations in DBI 
inflation. 
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